Multimodal single-cell analyses reveal mechanisms of perianal fistula in diverse patients with Crohn's disease.

BACKGROUND: Crohn's disease complicated by perianal fistulae is more prevalent and severe in patients of African ancestry. METHODS: We profiled single cells from diverse patients with Crohn's disease with perianal fistula from colorectal mucosa and fistulous tracts. Immunofluorescence was performed to validate predicted cell-cell interactions. Unstimulated monocytes were chronically cultured in diverse cohorts. A subset was analyzed by single-nucleus RNA + ATAC sequencing. FINDINGS: Fistulous tract cells from complete proctectomies demonstrated enrichment of myeloid cells compared to paired rectal tissues. Ligand-receptor analysis highlights myeloid-stromal cross-talk and cellular senescence, with cellular co-localization validated by immunofluorescence. Chitinase-3 like-protein-1 (CHI3L1) is a top upregulated gene in stromal cells from fistulae expressing both destructive and fibrotic gene signatures. Monocyte cultures from patients of African ancestry and controls demonstrated differences in CHI3L1 and oncostatin M (OSM) expression upon differentiation compared to individuals of European ancestry. Activating protein-1 footprints are present in ATAC-seq peaks in stress response genes, including CHI3L1 and OSM; genome-wide chromatin accessibility including JUN footprints was observed, consistent with reported mechanisms of inflammatory memory. Regulon analyses confirm known cell-specific transcription factor regulation and implicate novel ones in fibroblast subsets. All pseudo-bulked clusters demonstrate enrichment of genetic loci, establishing multicellular contributions. In the most significant African American Crohn's genetic locus, upstream of prostaglandin E receptor 4, lymphoid-predominant ATAC-seq peaks were observed, with predicted RORC footprints. CONCLUSIONS: Population differences in myeloid-stromal cross-talk implicate fibrotic and destructive fibroblasts, senescence, epigenetic memory, and cell-specific enhancers in perianal fistula pathogenesis. The transcriptomic and epigenetic data provided here may guide optimization of promising mesenchymal stem cell therapies for perianal fistula. FUNDING: This work was supported by grants U01DK062422, U24DK062429, and R01DK123758.

Gut-associated lymphoid tissue attrition associates with response to anti-α4ß7 therapy in ulcerative colitis.

Vedolizumab (VDZ) is a first-line treatment in ulcerative colitis (UC) that targets the α4ß7- mucosal vascular addressin cell adhesion molecule 1 (MAdCAM-1) axis. To determine the mechanisms of action of VDZ, we examined five distinct cohorts of patients with UC. A decrease in naïve B and T cells in the intestines and gut-homing (ß7+) plasmablasts in circulation of VDZ-treated patients suggested that VDZ targets gut-associated lymphoid tissue (GALT). Anti-α4ß7 blockade in wild-type and photoconvertible (KikGR) mice confirmed a loss of GALT size and cellularity because of impaired cellular entry. In VDZ-treated patients with UC, treatment responders demonstrated reduced intestinal lymphoid aggregate size and follicle organization and a reduction of ß7+IgG+ plasmablasts in circulation, as well as IgG+ plasma cells and FcγR-dependent signaling in the intestine. GALT targeting represents a previously unappreciated mechanism of action of α4ß7-targeted therapies, with major implications for this therapeutic paradigm in UC.

Gene coexpression networks reveal a broad role for lncRNAs in inflammatory bowel disease.

The role of long noncoding RNAs (lncRNAs) in disease is incompletely understood, but their regulation of inflammation is increasingly appreciated. We addressed the extent of lncRNA involvement in inflammatory bowel disease (IBD) using biopsy-derived RNA-sequencing data from a large cohort of deeply phenotyped patients with IBD. Weighted gene correlation network analysis revealed gene modules of lncRNAs coexpressed with protein-coding genes enriched for biological pathways, correlated with epithelial and immune cell signatures, or correlated with distal colon expression. Correlation of modules with clinical features uncovered a module correlated with disease severity, with an enriched interferon response signature containing the hub lncRNA IRF1-AS1. Connecting genes to IBD-associated single nucleotide polymorphisms (SNPs) revealed an enrichment of SNP-adjacent lncRNAs in biologically relevant modules. Ulcerative colitis-specific SNPs were enriched in distal colon-related modules, suggesting that disease-specific mechanisms may result from altered lncRNA expression. The function of the IBD-associated SNP-adjacent lncRNA IRF1-AS1 was explored in human myeloid cells, and our results suggested IRF1-AS1 promoted optimal production of TNF-α, IL-6, and IL-23. A CRISPR/Cas9-mediated activation screen in THP-1 cells revealed several lncRNAs that modulated LPS-induced TNF-α responses. Overall, this study uncovered the expression patterns of lncRNAs in IBD that identify functional, disease-relevant lncRNAs.

Integration of transcriptomes of senescent cell models with multi-tissue patient samples reveals reduced COL6A3 as an inducer of senescence.

Senescent cells are a major contributor to age-dependent cardiovascular tissue dysfunction, but knowledge of their in vivo cell markers and tissue context is lacking. To reveal tissue-relevant senescence biology, we integrate the transcriptomes of 10 experimental senescence cell models with a 224 multi-tissue gene co-expression network based on RNA-seq data of seven tissues biopsies from ∼600 coronary artery disease (CAD) patients. We identify 56 senescence-associated modules, many enriched in CAD GWAS genes and correlated with cardiometabolic traits-which supports universality of senescence gene programs across tissues and in CAD. Cross-tissue network analyses reveal 86 candidate senescence-associated secretory phenotype (SASP) factors, including COL6A3. Experimental knockdown of COL6A3 induces transcriptional changes that overlap the majority of the experimental senescence models, with cell-cycle arrest linked to modulation of DREAM complex-targeted genes. We provide a transcriptomic resource for cellular senescence and identify candidate biomarkers, SASP factors, and potential drivers of senescence in human tissues.

Myeloid cell influx into the colonic epithelium is associated with disease severity and non-response to anti-Tumor Necrosis Factor Therapy in patients with Ulcerative Colitis.

Ulcerative colitis (UC) is an idiopathic chronic inflammatory disease of the colon with sharply rising global prevalence. Dysfunctional epithelial compartment (EC) dynamics are implicated in UC pathogenesis although EC-specific studies are sparse. Applying orthogonal high-dimensional EC profiling to a Primary Cohort (PC; n=222), we detail major epithelial and immune cell perturbations in active UC. Prominently, reduced frequencies of mature BEST4+OTOP2+ absorptive and BEST2+WFDC2+ secretory epithelial enterocytes were associated with the replacement of homeostatic, resident TRDC+KLRD1+HOPX+ γδ+ T cells with RORA+CCL20+S100A4+ TH17 cells and the influx of inflammatory myeloid cells. The EC transcriptome (exemplified by S100A8, HIF1A, TREM1, CXCR1) correlated with clinical, endoscopic, and histological severity of UC in an independent validation cohort (n=649). Furthermore, therapeutic relevance of the observed cellular and transcriptomic changes was investigated in 3 additional published UC cohorts (n=23, 48 and 204 respectively) to reveal that non-response to anti-Tumor Necrosis Factor (anti-TNF) therapy was associated with EC related myeloid cell perturbations. Altogether, these data provide high resolution mapping of the EC to facilitate therapeutic decision-making and personalization of therapy in patients with UC.

Gut-associated lymphoid tissue attrition associates with response to anti-α4ß7 therapy in ulcerative colitis.

Targeting the α4ß7-MAdCAM-1 axis with vedolizumab (VDZ) is a front-line therapeutic paradigm in ulcerative colitis (UC). However, mechanism(s) of action (MOA) of VDZ remain relatively undefined. Here, we examined three distinct cohorts of patients with UC (n=83, n=60, and n=21), to determine the effect of VDZ on the mucosal and peripheral immune system. Transcriptomic studies with protein level validation were used to study drug MOA using conventional and transgenic murine models. We found a significant decrease in colonic and ileal naïve B and T cells and circulating gut-homing plasmablasts (ß7+) in VDZ-treated patients, pointing to gut-associated lymphoid tissue (GALT) targeting by VDZ. Murine Peyer's patches (PP) demonstrated a significant loss cellularity associated with reduction in follicular B cells, including a unique population of epithelium-associated B cells, following anti-α4ß7 antibody (mAb) administration. Photoconvertible (KikGR) mice unequivocally demonstrated impaired cellular entry into PPs in anti-α4ß7 mAb treated mice. In VDZ-treated, but not anti-tumor necrosis factor-treated UC patients, lymphoid aggregate size was significantly reduced in treatment responders compared to non-responders, with an independent validation cohort further confirming these data. GALT targeting represents a novel MOA of α4ß7-targeted therapies, with major implications for this therapeutic paradigm in UC, and for the development of new therapeutic strategies.

A temporal classifier predicts histopathology state and parses acute-chronic phasing in inflammatory bowel disease patients.

Previous studies have conducted time course characterization of murine colitis models through transcriptional profiling of differential expression. We characterize the transcriptional landscape of acute and chronic models of dextran sodium sulfate (DSS) and adoptive transfer (AT) colitis to derive temporal gene expression and splicing signatures in blood and colonic tissue in order to capture dynamics of colitis remission and relapse. We identify sub networks of patient-derived causal networks that are enriched in these temporal signatures to distinguish acute and chronic disease components within the broader molecular landscape of IBD. The interaction between the DSS phenotype and chronological time-point naturally defines parsimonious temporal gene expression and splicing signatures associated with acute and chronic phases disease (as opposed to ordinary time-specific differential expression/splicing). We show these expression and splicing signatures are largely orthogonal, i.e. affect different genetic bodies, and that using machine learning, signatures are predictive of histopathological measures from both blood and intestinal data in murine colitis models as well as an independent cohort of IBD patients. Through access to longitudinal multi-scale profiling from disease tissue in IBD patient cohorts, we can apply this machine learning pipeline to generation of direct patient temporal multimodal regulatory signatures for prediction of histopathological outcomes.

Biopsy and blood-based molecular biomarker of inflammation in IBD.

OBJECTIVE: IBD therapies and treatments are evolving to deeper levels of remission. Molecular measures of disease may augment current endpoints including the potential for less invasive assessments. DESIGN: Transcriptome analysis on 712 endoscopically defined inflamed (Inf) and 1778 non-inflamed (Non-Inf) intestinal biopsies (n=498 Crohn's disease, n=421 UC and 243 controls) in the Mount Sinai Crohn's and Colitis Registry were used to identify genes differentially expressed between Inf and Non-Inf biopsies and to generate a molecular inflammation score (bMIS) via gene set variance analysis. A circulating MIS (cirMIS) score, reflecting intestinal molecular inflammation, was generated using blood transcriptome data. bMIS/cirMIS was validated as indicators of intestinal inflammation in four independent IBD cohorts. RESULTS: bMIS/cirMIS was strongly associated with clinical, endoscopic and histological disease activity indices. Patients with the same histologic score of inflammation had variable bMIS scores, indicating that bMIS describes a deeper range of inflammation. In available clinical trial data sets, both scores were responsive to IBD treatment. Despite similar baseline endoscopic and histologic activity, UC patients with lower baseline bMIS levels were more likely treatment responders compared with those with higher levels. Finally, among patients with UC in endoscopic and histologic remission, those with lower bMIS levels were less likely to have a disease flare over time. CONCLUSION: Transcriptionally based scores provide an alternative objective and deeper quantification of intestinal inflammation, which could augment current clinical assessments used for disease monitoring and have potential for predicting therapeutic response and patients at higher risk of disease flares.

NOX1 is essential for TNFα-induced intestinal epithelial ROS secretion and inhibits M cell signatures.

OBJECTIVE: Loss-of-function mutations in genes generating reactive oxygen species (ROS), such as NOX1, are associated with IBD. Mechanisms whereby loss of ROS drive IBD are incompletely defined. DESIGN: ROS measurements and single-cell transcriptomics were performed on colonoids stratified by NOX1 genotype and TNFα stimulation. Clustering of epithelial cells from human UC (inflamed and uninflamed) scRNASeq was performed. Validation of M cell induction was performed by immunohistochemistry using UEA1 (ulex europaeus agglutin-1 lectin) and in vivo with DSS injury. RESULTS: TNFα induces ROS production more in NOX1-WT versus NOX1-deficient murine colonoids under a range of Wnt-mediated and Notch-mediated conditions. scRNASeq from inflamed and uninflamed human colitis versus TNFα stimulated, in vitro colonoids defines substantially shared, induced transcription factors; NOX1-deficient colonoids express substantially lower levels of STAT3 (signal transducer and activator of transcription 3), CEBPD (CCAAT enhancer-binding protein delta), DNMT1 (DNA methyltransferase) and HIF1A (hypoxia-inducible factor) baseline. Subclustering unexpectedly showed marked TNFα-mediated induction of M cells (sentinel cells overlying lymphoid aggregates) in NOX1-deficient colonoids. M cell induction by UEA1 staining is rescued with H2O2 and paraquat, defining extra- and intracellular ROS roles in maintenance of LGR5+ stem cells. DSS injury demonstrated GP2 (glycoprotein-2), basal lymphoplasmacytosis and UEA1 induction in NOX1-deficiency. Principal components analyses of M cell genes and decreased DNMT1 RNA velocity correlate with UC inflammation. CONCLUSIONS: NOX1 deficiency plus TNFα stimulation contribute to colitis through dysregulation of the stem cell niche and altered cell differentiation, enhancing basal lymphoplasmacytosis. Our findings prioritise ROS modulation for future therapies.

T3 and glucose increase expression of phosphoenolpyruvate carboxykinase (PCK1) leading to increased ß-cell proliferation.

OBJECTIVES: Thyroid hormone (T3) and high glucose concentrations are critical components of ß-cell maturation and function. In the present study, we asked whether T3 and glucose signaling pathways coordinately regulate transcription of genes important for ß-cell function and proliferation. METHODS: RNA-seq analysis was performed on cadaveric human islets from five different donors in response to low and high glucose concentrations and in the presence or absence of T3. Gene expression was also studies in sorted human ß-cells, mouse islets and Ins-1 cells by RT-qPCR. Silencing of the thyroid hormone receptors (THR) was conducted using lentiviruses. Proliferation was assessed by ki67 immunostaining in primary human/mouse islets. Chromatin immunoprecipitation and proximity ligation assay were preformed to validate interactions of ChREBP and THR. RESULTS: We found glucose-mediated expression of carbohydrate response element binding protein alpha and beta (ChREBPα and ChREBPß) mRNAs and their target genes are highly dependent on T3 concentrations in rodent and human ß-cells. In ß-cells, T3 and glucose coordinately regulate the expression of ChREBPß and PCK1 (phosphoenolpyruvate carboxykinase-1) among other important genes for ß-cell maturation. Additionally, we show the thyroid hormone receptor (THR) and ChREBP interact, and their relative response elements are located near to each other on mutually responsive genes. In FACS-sorted adult human ß-cells, we found that high concentrations of glucose and T3 induced the expression of PCK1. Next, we show that overexpression of Pck1 together with dimethyl malate (DMM), a substrate precursor, significantly increased ß-cell proliferation in human islets. Finally, using a Cre-Lox approach, we demonstrated that ChREBPß contributes to Pck1-dependent ß-cell proliferation in mouse ß-cells. CONCLUSIONS: We conclude that T3 and glucose act together to regulate ChREBPß, leading to increased expression and activity of Pck1, and ultimately increased ß-cell proliferation.

A mechanistic framework for cardiometabolic and coronary artery diseases.

Coronary atherosclerosis results from the delicate interplay of genetic and exogenous risk factors, principally taking place in metabolic organs and the arterial wall. Here we show that 224 gene-regulatory coexpression networks (GRNs) identified by integrating genetic and clinical data from patients with (n = 600) and without (n = 250) coronary artery disease (CAD) with RNA-seq data from seven disease-relevant tissues in the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) study largely capture this delicate interplay, explaining >54% of CAD heritability. Within 89 cross-tissue GRNs associated with clinical severity of CAD, 374 endocrine factors facilitated inter-organ interactions, primarily along an axis from adipose tissue to the liver (n = 152). This axis was independently replicated in genetically diverse mouse strains and by injection of recombinant forms of adipose endocrine factors (EPDR1, FCN2, FSTL3 and LBP) that markedly altered blood lipid and glucose levels in mice. Altogether, the STARNET database and the associated GRN browser (http://starnet.mssm.edu) provide a multiorgan framework for exploration of the molecular interplay between cardiometabolic disorders and CAD.

Disrupting the DREAM complex enables proliferation of adult human pancreatic ß cells.

Resistance to regeneration of insulin-producing pancreatic ß cells is a fundamental challenge for type 1 and type 2 diabetes. Recently, small molecule inhibitors of the kinase DYRK1A have proven effective in inducing adult human ß cells to proliferate, but their detailed mechanism of action is incompletely understood. We interrogated our human insulinoma and ß cell transcriptomic databases seeking to understand why ß cells in insulinomas proliferate, while normal ß cells do not. This search reveals the DREAM complex as a central regulator of quiescence in human ß cells. The DREAM complex consists of a module of transcriptionally repressive proteins that assemble in response to DYRK1A kinase activity, thereby inducing and maintaining cellular quiescence. In the absence of DYRK1A, DREAM subunits reassemble into the pro-proliferative MMB complex. Here, we demonstrate that small molecule DYRK1A inhibitors induce human ß cells to replicate by converting the repressive DREAM complex to its pro-proliferative MMB conformation.

A mitochondrial long-chain fatty acid oxidation defect leads to transfer RNA uncharging and activation of the integrated stress response in the mouse heart.

AIMS: Cardiomyopathy and arrhythmias can be severe presentations in patients with inherited defects of mitochondrial long-chain fatty acid ß-oxidation (FAO). The pathophysiological mechanisms that underlie these cardiac abnormalities remain largely unknown. We investigated the molecular adaptations to a FAO deficiency in the heart using the long-chain acyl-CoA dehydrogenase (LCAD) knockout (KO) mouse model. METHODS AND RESULTS: We observed enrichment of amino acid metabolic pathways and of ATF4 target genes among the upregulated genes in the LCAD KO heart transcriptome. We also found a prominent activation of the eIF2α/ATF4 axis at the protein level that was independent of the feeding status, in addition to a reduction of cardiac protein synthesis during a short period of food withdrawal. These findings are consistent with an activation of the integrated stress response (ISR) in the LCAD KO mouse heart. Notably, charging of several transfer RNAs (tRNAs), such as tRNAGln was decreased in LCAD KO hearts, reflecting a reduced availability of cardiac amino acids, in particular, glutamine. We replicated the activation of the ISR in the hearts of mice with muscle-specific deletion of carnitine palmitoyltransferase 2. CONCLUSIONS: Our results show that perturbations in amino acid metabolism caused by long-chain FAO deficiency impact cardiac metabolic signalling, in particular the ISR. These results may serve as a foundation for investigating the role of the ISR in the cardiac pathology associated with long-chain FAO defects.Translational Perspective: The heart relies mainly on mitochondrial fatty acid ß-oxidation (FAO) for its high energy requirements. The heart disease observed in patients with a genetic defect in this pathway highlights the importance of FAO for cardiac health. We show that the consequences of a FAO defect extend beyond cardiac energy homeostasis and include amino acid metabolism and associated signalling pathways such as the integrated stress response.

Ulcerative colitis is characterized by a plasmablast-skewed humoral response associated with disease activity.

B cells, which are critical for intestinal homeostasis, remain understudied in ulcerative colitis (UC). In this study, we recruited three cohorts of patients with UC (primary cohort, n = 145; validation cohort 1, n = 664; and validation cohort 2, n = 143) to comprehensively define the landscape of B cells during UC-associated intestinal inflammation. Using single-cell RNA sequencing, single-cell IgH gene sequencing and protein-level validation, we mapped the compositional, transcriptional and clonotypic landscape of mucosal and circulating B cells. We found major perturbations within the mucosal B cell compartment, including an expansion of naive B cells and IgG+ plasma cells with curtailed diversity and maturation. Furthermore, we isolated an auto-reactive plasma cell clone targeting integrin αvß6 from inflamed UC intestines. We also identified a subset of intestinal CXCL13-expressing TFH-like T peripheral helper cells that were associated with the pathogenic B cell response. Finally, across all three cohorts, we confirmed that changes in intestinal humoral immunity are reflected in circulation by the expansion of gut-homing plasmablasts that correlates with disease activity and predicts disease complications. Our data demonstrate a highly dysregulated B cell response in UC and highlight a potential role of B cells in disease pathogenesis.

An Integrated Taxonomy for Monogenic Inflammatory Bowel Disease.

BACKGROUND & AIMS: Monogenic forms of inflammatory bowel disease (IBD) illustrate the essential roles of individual genes in pathways and networks safeguarding immune tolerance and gut homeostasis. METHODS: To build a taxonomy model, we assessed 165 disorders. Genes were prioritized based on penetrance of IBD and disease phenotypes were integrated with multi-omics datasets. Monogenic IBD genes were classified by (1) overlapping syndromic features, (2) response to hematopoietic stem cell transplantation, (3) bulk RNA-sequencing of 32 tissues, (4) single-cell RNA-sequencing of >50 cell subsets from the intestine of healthy individuals and patients with IBD (pediatric and adult), and (5) proteomes of 43 immune subsets. The model was validated by addition of newly identified monogenic IBD defects. As a proof-of-concept, we explore the intersection between immunometabolism and antimicrobial activity for a group of disorders (G6PC3/SLC37A4). RESULTS: Our quantitative integrated taxonomy defines the cellular landscape of monogenic IBD gene expression across 102 genes with high and moderate penetrance (81 in the model set and 21 genes in the validation set). We illustrate distinct cellular networks, highlight expression profiles across understudied cell types (e.g., CD8+ T cells, neutrophils, epithelial subsets, and endothelial cells) and define genotype-phenotype associations (perianal disease and defective antimicrobial activity). We illustrate processes and pathways shared across cellular compartments and phenotypic groups and highlight cellular immunometabolism with mammalian target of rapamycin activation as one of the converging pathways. There is an overlap of genes and enriched cell-specific expression between monogenic and polygenic IBD. CONCLUSION: Our taxonomy integrates genetic, clinical and multi-omic data; providing a basis for genomic diagnostics and testable hypotheses for disease functions and treatment responses.

Integrative Analysis of the Inflammatory Bowel Disease Serum Metabolome Improves Our Understanding of Genetic Etiology and Points to Novel Putative Therapeutic Targets.

BACKGROUND & AIMS: Polygenic and environmental factors are underlying causes of inflammatory bowel disease (IBD). We hypothesized that integration of the genetic loci controlling a metabolite's abundance, with known IBD genetic susceptibility loci, may help resolve metabolic drivers of IBD. METHODS: We measured the levels of 1300 metabolites in the serum of 484 patients with ulcerative colitis (UC) and 464 patients with Crohn's disease (CD) and 365 controls. Differential metabolite abundance was determined for disease status, subtype, clinical and endoscopic disease activity, as well as IBD phenotype including disease behavior, location, and extent. To inform on the genetic basis underlying metabolic diversity, we integrated metabolite and genomic data. Genetic colocalization and Mendelian randomization analyses were performed using known IBD risk loci to explore whether any metabolite was causally associated with IBD. RESULTS: We found 173 genetically controlled metabolites (metabolite quantitative trait loci, 9 novel) within 63 non-overlapping loci (7 novel). Furthermore, several metabolites significantly associated with IBD disease status and activity as defined using clinical and endoscopic indexes. This constitutes a resource for biomarker discovery and IBD biology insights. Using this resource, we show that a novel metabolite quantitative trait locus for serum butyrate levels containing ACADS was not supported as causal for IBD; replicate the association of serum omega-6 containing lipids with the fatty acid desaturase 1/2 locus and identify these metabolites as causal for CD through Mendelian randomization; and validate a novel association of serum plasmalogen and TMEM229B, which was predicted as causal for CD. CONCLUSIONS: An exploratory analysis combining genetics and unbiased serum metabolome surveys can reveal novel biomarkers of disease activity and potential mediators of pathology in IBD.

Peroxisomal L-bifunctional protein (EHHADH) deficiency causes male-specific kidney hypertrophy and proximal tubular injury in mice.

BACKGROUND: Proximal tubular (PT) cells are enriched in mitochondria and peroxisomes. Whereas mitochondrial fatty acid oxidation (FAO) plays an important role in kidney function by supporting the high-energy requirements of PT cells, the role of peroxisomal metabolism remains largely unknown. EHHADH, also known as L-bifunctional protein, catalyzes the second and third step of peroxisomal FAO. METHODS: We studied kidneys of WT and Ehhadh KO mice on a C57BL/6N background using histology, immunohistochemistry, immunofluorescence, immunoblot, RNA-sequencing, and metabolomics. To assess the role of androgens in the kidney phenotype of Ehhadh KO mice, mice underwent orchiectomy. RESULTS: We observed male-specific kidney hypertrophy and glomerular filtration rate reduction in adult Ehhadh KO mice. Transcriptome analysis unveiled a gene expression signature similar to PT injury in acute kidney injury mouse models. This was further illustrated by the presence of KIM-1 (kidney injury molecule-1), SOX-9, and Ki67-positive cells in the PT of male Ehhadh KO kidneys. Male Ehhadh KO kidneys had metabolite changes consistent with peroxisomal dysfunction as well as an elevation in glycosphingolipid levels. Orchiectomy of Ehhadh KO mice decreased the number of KIM-1 positive cells to WT levels. We revealed a pronounced sexual dimorphism in the expression of peroxisomal FAO proteins in mouse kidney, underlining a role of androgens in the kidney phenotype of Ehhadh KO mice. CONCLUSIONS: Our data highlight the importance of EHHADH and peroxisomal metabolism in male kidney physiology and reveal peroxisomal FAO as a sexual dimorphic metabolic pathway in mouse kidneys.

Molecular Characterization of Limited Ulcerative Colitis Reveals Novel Biology and Predictors of Disease Extension.

BACKGROUND AND AIMS: Disease extent varies in ulcerative colitis (UC) from proctitis to left-sided colitis to pancolitis and is a major prognostic factor. When the extent of UC is limited there is often a sharp demarcation between macroscopically involved and uninvolved areas and what defines this or subsequent extension is unknown. We characterized the demarcation site molecularly and determined genes associated with subsequent disease extension. METHODS: We performed RNA sequence analysis of biopsy specimens from UC patients with endoscopically and histologically confirmed limited disease, of which a subset later extended. Biopsy specimens were obtained from the endoscopically inflamed upper (proximal) limit of disease, immediately adjacent to the uninvolved colon, as well as at more proximal, endoscopically uninflamed colonic segments. RESULTS: Differentially expressed genes were identified in the endoscopically inflamed biopsy specimens taken at each patient's most proximal diseased site relative to healthy controls. Expression of these genes in the more proximal biopsy specimens transitioned back to control levels abruptly or gradually, the latter pattern supporting the concept that disease exists beyond the endoscopic disease demarcation site. The gradually transitioning genes were associated with inflammation, angiogenesis, glucuronidation, and homeodomain pathways. A subset of these genes in inflamed biopsy specimens was found to predict disease extension better than clinical features and were responsive to biologic therapies. Network analysis revealed critical roles for interferon signaling in UC inflammation and poly(ADP-ribose) polymerase 14 (PARP14) was a predicted key driver gene of extension. Higher PARP14 protein levels were found in inflamed biopsy specimens of patients with limited UC that subsequently extended. CONCLUSION: Molecular predictors of disease extension reveal novel strategies for disease prognostication and potential therapeutic targeting.

Stratification of risk of progression to colectomy in ulcerative colitis via measured and predicted gene expression.

An important goal of clinical genomics is to be able to estimate the risk of adverse disease outcomes. Between 5% and 10% of individuals with ulcerative colitis (UC) require colectomy within 5 years of diagnosis, but polygenic risk scores (PRSs) utilizing findings from genome-wide association studies (GWASs) are unable to provide meaningful prediction of this adverse status. By contrast, in Crohn disease, gene expression profiling of GWAS-significant genes does provide some stratification of risk of progression to complicated disease in the form of a transcriptional risk score (TRS). Here, we demonstrate that a measured TRS based on bulk rectal gene expression in the PROTECT inception cohort study has a positive predictive value approaching 50% for colectomy. Single-cell profiling demonstrates that the genes are active in multiple diverse cell types from both the epithelial and immune compartments. Expression quantitative trait locus (QTL) analysis identifies genes with differential effects at baseline and week 52 follow-up, but for the most part, differential expression associated with colectomy risk is independent of local genetic regulation. Nevertheless, a predicted polygenic transcriptional risk score (PPTRS) derived by summation of transcriptome-wide association study (TWAS) effects identifies UC-affected individuals at 5-fold elevated risk of colectomy with data from the UK Biobank population cohort studies, independently replicated in an NIDDK-IBDGC dataset. Prediction of gene expression from relatively small transcriptome datasets can thus be used in conjunction with TWASs for stratification of risk of disease complications.

Murine deficiency of peroxisomal L-bifunctional protein (EHHADH) causes medium-chain 3-hydroxydicarboxylic aciduria and perturbs hepatic cholesterol homeostasis.

Peroxisomes play an essential role in the ß-oxidation of dicarboxylic acids (DCAs), which are metabolites formed upon ω-oxidation of fatty acids. Genetic evidence linking transporters and enzymes to specific DCA ß-oxidation steps is generally lacking. Moreover, the physiological functions of DCA metabolism remain largely unknown. In this study, we aimed to characterize the DCA ß-oxidation pathway in human cells, and to evaluate the biological role of DCA metabolism using mice deficient in the peroxisomal L-bifunctional protein (Ehhadh KO mice). In vitro experiments using HEK-293 KO cell lines demonstrate that ABCD3 and ACOX1 are essential in DCA ß-oxidation, whereas both the bifunctional proteins (EHHADH and HSD17B4) and the thiolases (ACAA1 and SCPx) have overlapping functions and their contribution may depend on expression level. We also show that medium-chain 3-hydroxydicarboxylic aciduria is a prominent feature of EHHADH deficiency in mice most notably upon inhibition of mitochondrial fatty acid oxidation. Using stable isotope tracing methodology, we confirmed that products of peroxisomal DCA ß-oxidation can be transported to mitochondria for further metabolism. Finally, we show that, in liver, Ehhadh KO mice have increased mRNA and protein expression of cholesterol biosynthesis enzymes with decreased (in females) or similar (in males) rate of cholesterol synthesis. We conclude that EHHADH plays an essential role in the metabolism of medium-chain DCAs and postulate that peroxisomal DCA ß-oxidation is a regulator of hepatic cholesterol biosynthesis.